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Foreground subtraction is the biggest challenge for future redshifted 21cm ob- 
servations to probe reionization. We use a short GMRT observation at 153 MHz to 
characterize the statistical properties of the background radiation across ~ 1° to 
sub-arcminutes angular scales, and across a frequency band of 5 MHz with 62.5 kHz 
resolution. The statistic we use is the visibility correlation function, or cquivalently 
the angular power spectrum Ci . We present the results obtained from using relatively 
unsophisticated, conventional data calibration procedures. We find that even fairly 
simple minded calibration allows one to estimate the visibility correlation function at 
a given frequency Vz(U, 0). From our observations we find that Vz(U,0) is consistent 
with foreground model predictions at all angular scales except the largest ones probed 
by our observations where the the model predictions are somewhat in excess. On the 
other hand the visibility correlation between different frequencies k(U,Av), seems to 
be much more sensitive to calibration errors. We find a rapid decline in k(U, Av), in 
contrast with the prediction of less than 1% variation across 2.5 MHz. In this case 
however, it seems likely that a substantial part of the discrepancy may be due to 
limitations of data reduction procedures. 

Key words: cosmology: observations, cosmology: diffuse radiation, methods: statis- 
tical 
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1 INTRODUCTION 

Observations of redshifted 21cm radiation from the large 
scale distribution of neutral hydrogen (HI) are perceived as 
one of the most promising future probes of the Universe 
at high redshifts (see Furlanetto ,Oh & Briggs 2006 for a 
recent review). Observational evidence from quasar absorp- 
tion spectra (Becker et al. 2001; Fan et al. 2002) and the 
CMBR (Spergel et al. 2007; Page et al. 2007) together im- 
ply that the HI was reionized over an extended period span- 
ning the redshift range 6 < z < 15 (for reviews see Barkana 
& Loeb 2001; Fan, Carilli & Keating 2006; Choudhury & 
Ferrara 2006). Determining how and when the Universe was 
reionized is one of the most important issues that will be 
addressed by future 21cm observations. The Giant Meter 
Wave Radio Telescope (GMRT l ; Swarup et al. 1991), cur- 
rently functioning at several frequency bands in the range 
150 to 1420 MHz is very well suited for carrying out initial 



investigations towards detecting the reionization HI signal. 
There are several upcoming low-frequency instruments such 
as LOFAR 2 , MWA 3 , 21CMA 4 and SKA 5 which are being 
built specifically with these observations in view. 

It is currently perceived that a statistical analysis of the 
fluctuations in the redshifted 21 cm signal holds the great- 
est potential for observing HI at high redshifts (Bharad- 
waj and Sethi 2001; Zaldarriaga, Furlanetto & Hernquist 
2004; Morales & Hewitt 2004; Bharadwaj and Ali 2005; 
Bharadwaj and Pandey 2005). Correlations among the vis- 
ibilities measured in radio-interferometric observations di- 
rectly probe the HI power spectrum at the epoch where the 
radiation originated. The reionization visibility signal at the 
GMRT is expected to be ~ 1 mjy and smaller (Bharadwaj 
and Ali 2005). This HI signal is present as a minute com- 
ponent of the background in all low frequency observations, 
and it is buried in foreground radiation from other astro- 
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physical sources whose contribution is 4 to 5 orders of mag- 
nitude larger. Extracting the HI signal from the foregrounds 
is a major challenge. 

Individual sources can be identified and removed from 
the image at a flux level which depends on the sensitivity. 
The contribution from the remaining discrete sources could 
be large enough to overwhelm the HI signal (Di Matteo et 
al. 2002) . The diffuse synchrotron emission from our Galaxy 
(Shaver et al. 1999) is another important component. Fore- 
ground sources include free-free emission from ionizing halos 
(Oh & Mack 2003), faint radio loud quasars (Di Matteo et 
al. 2002) and synchrotron emission from low redshift galaxy 
clusters (DiMatteo. et al. 2004). 

The foregrounds are expected to have a continuum spec- 
tra, and the contribution at two different frequencies sepa- 
rated by Af ~ 1 MHz are expected to be highly correlated. 
The HI signal is expected to be uncorrelated at such a fre- 
quency separation and this holds the promise of allowing us 
to separate the signal from the foregrounds. A possible ap- 
proach is to subtract a best fit continuum spectra for each 
line of sight (Wang et al. 2006) and then use the residuals to 
determine the HI power spectrum. An alternate approach is 
to first determine the statistical properties of the total radi- 
ation and then subtract out the smooth Ai; dependent part 
to extract the HI signal (Zaldarriaga, Furlanetto & Hern- 
quist 2004). The issue of foreground removal has also been 
studied by Morales et al. (2006) and Mcquinn et al. (2006). 

It is crucial to accurately characterize the foregrounds 
in order to be able to detect the HI signal in future observa- 
tions. In this paper we used GMRT observations to charac- 
terize the foregrounds at 153 MHz which corresponds to an 
HI signal from z = 8.3. To the best of our knowledge this is 
the first attempt to directly characterize the foregrounds at 
angular scales (~ 1° to sub-arcminute) and frequency cover- 
age (6 MHz with 62.5 kHz resolution) relevant for detecting 
the reionization HI signal. 

We next present a brief outline of the paper. In Sec- 
tion 2 we describe the observations and data reduction 
while in Section 3 we present "visibility-correlations" which 
we use to quantify the statistical properties of our radio- 
interferometric data. Section 4 presents the predictions of 
existing foreground models, and in Section 5 we present our 
results and discuss their implications. 



2 GMRT OBSERVATIONS AND DATA 
REDUCTION 

The GMRT has a hybrid configuration (Swarup et al. 1991) 
where 14 of the 30 antennas are randomly distributed in a 
Central Square ~ 1.1km x 1.1km in extent. These provide 
the uv coverage at small baselines. Here baseline refers to the 
antenna separation, and we use the two dimensional vector 
U to denote the component perpendicular to the direction 
of observation. Note that U has Cartesian components (u, v) 
and is dimensionless being in units of the observing wave- 
length. The shortest baseline at the GMRT is 100 m which 
comes down to around 60 m with projection effects. The rest 
of the antennas in the GMRT lie along three arms in an ap- 
proximately 'Y' configuration. These provide uv coverage at 
long baselines (the longest baseline is 26 km) . The diameter 
of each GMRT antenna is 45m. The hybrid configuration 
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Figure 1. This shows the uv coverage of the GMRT data that we 
have analyzed. Here (u, v) are the antenna separations in wave- 
length units at the observing frequency 153 MHz. 



gives reasonably good sensitivity for both compact and ex- 
tended sources. Figure 1 shows the uv coverage of our GMRT 
observations. 

On 15 th June, 2005 we observed a field centered on Up- 
silon Andromedae (which is an extra-solar planetary system 
system at a 2 ooo = l h 36 m 48", <5 20 oo = 41°24'23") for a total 
of 14 hours (including calibration). No emission that could 
be attributed to the planet was detected in our observations. 
The galactic co-ordinates are I = 132.00°, b = 20.67°. From 
the 408 MHz Haslam et al. (1982) map the sky temperature 
at this location is ~ 30 K (at 408 MHz), and there is no 
structure visible at the angular resolution of the map. 

The observational set up used a total of 128 frequency 
channels spanning 8 MHz centered at 153 MHz. Each fre- 
quency channel is 62.5kHz wide. A 6 MHz wide band-pass 
filter was introduced in the IF stage to exclude known strong 
Radio Frequency Interference (RFI), hence only 3/4 of the 
central channels contain astronomical signals. The integra- 
tion time was 16 seconds, and visibilities were recorded for 
two orthogonal circular polarizations. The visibility data 
were analyzed using the Astronomical Image Processing 
Software (AIPS). The calibrator source 3C48 was used for 
flux, phase and bandpass calibration. The calibrator was 
observed every half hour so as to correct for temporal vari- 
ations in the system gain. Standard AIPS tasks were used 
to flag all data that could be visually identified as being 
bad. We then made a high resolution image of the source 
using only a single channel (channel 35). The synthesized 
beam has a FWHM of 29" x 25" and the rms. noise in the 
CLEANed image is 9.5mJy/Beam. All sources with flux 
density more than 30mJy were fitted with clean compo- 
nents (CC), these components were merged and the visi- 
bilities corresponding to components with flux more than 
8.6 mjy were subtracted from the multi-channel uv data us- 
ing UVSUB. The value 8.6 mjy was chosen because we find 
predominantly positive clean components above this flux 
level whereas positive and negative components are equally 
abundant below this. The resultant uv data is now expected 
to be dominated by noise and residual RFI, since the major- 
ity of the point sources have been removed. Visually inspect- 
ing the data using the AIPS tasks VPLOT and UVHGM, 
we decided to clip the data at 12 Jy whereby visibilities with 
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amplitude greater than 12 Jy were discarded. The clipping 
amplitude is in principle crucial since one would like to en- 
sure that all baselines with RFI contributions have been dis- 
carded, without throwing away any good baselines. In prac- 
tice we found that the exact clipping value does not make a 
substantial difference in our subsequent analysis. After this 
we added back the visibilities corresponding to all the CC 
components that we had subtracted. To first order, one could 
expected that at this stage all strong RFI has been removed. 

The large field of view (6Whm = 3.8°) of the GMRT at 
150 MHz lead to considerable errors if the non-planar nature 
of the GMRT antenna distribution is not taken into account. 
We use the three dimensional (3D) imaging feature (e.g., 
Perley 1999) in the MPS task IMAGR in which the entire 
field of view is divided into multiple subfields (facets) each of 
which is imaged separately. Here a 4° x 4° field of view was 
imaged using 139 facets. We first collapsed 10 adjacent chan- 
nels (channels 30 to 39) to make a single channel which was 
used to make a CLEANed image. This channel's frequency 
width 0.625 MHz (< 0.7 MHz) which is sufficiently small 
so as to avoid bandwidth smearing. The synthesized beam 
has FWHM ~ 20" and the cleaned image has rms noise 
4.6mJy/Beam. The presence of a large number of sources 
in the field allows us to do self calibration loops to improve 
the image quality. The data went through 4 rounds of phase 
self calibration and a 5 th round where self calibration was 
done for both amplitude and phase. The time interval for the 
gain correction was chosen as 5, 5, 2, 2 and 2 minutes for the 
successive self calibration loops. The rms. noise in the final 
cleaned image was 3.1 mJy/Beam and the image quality had 
improved considerably. The final gain table was applied to 
all 128 frequency channels. Channels 21 to 100 of this data 
were then collapsed into 8 channels, each containing 10 of the 
original frequency channels. We use these to make a contin- 
uum image of the entire field. Some more data was flagged at 
this stage, and we then applied a final phase self calibration 
loop. This calibrated data was used to make the final cleaned 
image which is shown in Figure 2. The synthesized beam has 
a FWHM of 28" x 23", and an off-source RMS. noise level 
of 1.6mJy/Beam. Note that several of the extended features 
like the one at 02000 = 01 h 41 m , <5 2 ooo = 40°24 are actu- 
ally imaging artifacts around the brightest point sources. 
The brightest sources are also found to be accompanied 
by a region of negative flux density, these are presumably 
the results if residual phase errors which were not corrected 
for in our self calibration process. The maximum and mini- 
mum flux density in the final image are 820 mjy /Beam and 
— 44 m Jy /Beam respectively. 

Recall that for this experiment, the sources visible in 
the final continuum image (Figure 2) are contaminants 
which have to be removed. Pixels with flux density above 
8mJy/Beam which were visually identified as sources and 
not imaging artifacts were fitted with clean components. 
The clean components were merged and the visibilities cor- 
responding to these clean components were subtracted from 
the original full frequency resolution uv data using the AIPS 
task UVSUB. It is expected that at this stage most of 
the genuine sources in Figure 2 have been removed from 
the data. Figure 3 shows the final image made from the 
residual visibility data after UVSUB. The maximum and 
minimum flux density in this image are 25mJy/Beam and 
— 45mJy/Beam respectively. The subsequent analysis was 



done using the visibility data. We have analyzed the data 
both before and after the sources were subtracted, and we 
shall refer to these as data I (Initial - before source sub- 
traction) and data R (residual - after source subtraction) 
respectively. 

The final data contains 295868 baselines, each of which 
has visibilities for 2 circular polarizations and 96 frequency 
channels, of which we have used only the first 80 channels 
for the subsequent analysis. The visibilities from the two po- 
larizations were combined for the subsequent analysis. The 
real and imaginary parts of the resulting visibilities have a 
mean value —0.56 mjy and 2.6 mjy respectively , and rms of 
2.93 Jy for both in data I. For data R the real and imag- 
inary parts of the visibilities have a mean — 6.0mJy and 
1.1 mjy respectively whereas the rms is 2.42 Jy for both. 

In the subsequent analysis it is often convenient to as- 
sume that the visibilities have a Gaussian distribution. Fig- 
ure 4 shows the distribution of the real part of the visi- 
bilities for data R. We find that a Gaussian gives a rea- 
sonably good fit to the data within 2a which contains the 
bulk of the data. The number counts predicted by the Gaus- 
sian falls much faster than the data at large visibility values 
I Re(V) |> 6Jy. Deviation from Gaussian statistics is ex- 
pected to mainly affect the error estimate on the visibility 
correlation. We expect this effect to be small, since only 
a small fraction of visibilities are discrepant. The imaginary 
part of of data R, and the real and imaginary parts of data 
I all show a similar behaviour. 



3 VISIBILITY CORRELATIONS 

The visibility V(\J,v) measured in a radio-interferometric 
observation is the sum of three different contributions 

V{U,v) = S(\J,v) + F{U,v) + N(\J,v) (1) 

the HI signal S(U, v), astrophysical foregrounds F(U, v) 
and system noise N(U,v). We treat all three of these con- 
tributions as uncorrelated random variables with zero mean. 
The statistical properties of the visibility can be quantified 
through the two visibility correlation (henceforth the visi- 
bility correlation) 

^(Ui,i/i;U 2 ,i/2) = <V(Ui,z/i)V*(U 2 ,!/ 2 )} (2) 
and 

V 2 = S 2 + F 2 + N 2 (3) 

where S 2 , F 2 and N 2 respectively refer to the signal, 
foreground and noise contributions to the visibility correla- 
tion. 

The contribution from the HI signal S 2 is expected to 
be ~ 10~ 7 Jy 2 or smaller at 150 MHz (Bharadwaj and Ali 
2005). This is negligible compared to the expected fore- 
grounds and noise contributions in our observations, and 
hence we ignore it in our further analysis. 

The foreground contribution F(\J, v) is the Fourier 
transform of the product of the foreground specific inten- 
sity distribution on the sky 1(9, v) and the primary beam 
pattern of the individual GMRT antenna A(9, v). As men- 
tioned earlier, this Fourier relation is strictly valid only if 
the field of view is small, and in this observation we expect 
considerable deviations at large baselines. As we are mainly 
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Figure 2. These shows our continuum image of bandwidth 5 MHz centered at 153 MHz. The 4° X 4° field was imaged using 139 facets 
which have been combined using the AIPS task FLATN. The rms noise is 1.6mJy/Beam. The left and right panels shows positive and 
negative 7 — a contours respectively. Note that many of the extended positive features and all negative features are imaging artifacts 
around the brightest sources. 
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Cont peak flux = -4.4664E-02 JY/BEAM 
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Figure 3. This is the same as the Figure 2 except that all the bright pixels > 8mJy/Bcam that were visually identified as being genuine 
sources and not artifacts have been fitted with clean components and removed from the visibility data from which this image was made. 
It is expected that most of the genuine sources have been removed from this data. 



interested in the visibility correlations at small baselines, 
and also because the analysis is considerably more compli- 
cated otherwise, we assume the Fourier relation to hold. We 
can then express F(XJ, v) as a convolution 

F(U,z/) = J J(U>)a(U-U',i/)fi 2 U'. (4) 

where 7(U, v) and a(U, v) are the Fourier transform of 
1(9, v) and A(8, v) respectively. Assuming that the region 
of sky under observation is small so that it can be treated 



as flat, we have 

(/(Ur^/CU^)) = ^(Ur-U,)^)^)^ 

x CWe/iOi,^) (5) 

where <5f,(Ui — U2) is the two dimensional Dirac Delta func- 
tion, (dB/dT) v = 2k B v 2 /c 2 is the conversion factor from 
brightness temperature to specific intensity and Ci{v\,v-2) 
is the multi-frequency angular power spectrum (MAPS; eg. 
Datta, Roy Choudhury, & Bharadwaj 2007) of the fore- 
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Figure 4. The distribution of visibilities after source subtraction (data R). The same plot is shown on a linear scale (left panel) and 
a log-linear scale (right panel). The data is plotted as a histogram, and a Gaussian with the corresponding mean and rms. (see text) is 
plotted as a solid line. The discrepancy at high amplitudes (> 6 Jy) is visible only in the right panel. 



ground brightness temperature distribution. Using this to 
calculate the foreground contribution to the visibility corre- 
lation we have 

F2(Ui,i/i;U2,i*2) = ydVa(Ui-U>i)a*(U 2 -U> 2 ) 

>< (§?)„(£)„ <*> 

The GMRT primary beam is well parametrized by a Gaus- 
sian A(9,v) = e - fl2/e o where 9 w 0.6 x 6> FW hm = 2.3°. 
There is a small variation in do (oc v 1 ) across the fre- 
quency band. Ignoring this v dependence have 5(U, v) = 
o(U) = tt8q exp[— 9 2 % 2 U 2 ]. The integral in eq. (6) has a very 
small value unless the terms 5(Ui — U ) and fi*(U2 — U ) 
have a considerable overlap ie. | Ui — U2 |< (tt9q)~ ■ This 
tells us that i*2(U, vi;U + AU, 1/2) has a significant value 
only if |AU| < (tt9o)~ 1 and is negligible otherwise. Further, 
j AU !<C U at the baselines of interest, and we may ap- 
proximate a*(U + AU - U') as a*(U - U') in eq. (6) and 
write 

F 2 (XJ,v;XJ + AXJ,v + Av) = (H) 2 J dV|5(U-U')| 2 

x C 2nU ,(v,v + Av). (7) 

where we have ignored the Av dependence of 9o and 

The explicit reference to AU can be dropped as it does 
not appear in the integral. We also assume that CWt/^i, ^2) 
is a slowly varying function of U as compared to a(U)| 2 
whereby |a(U - U')| 2 w (it9 2 ] /2)S 2 d (U - U') which gives 

F 2 (U,Az,) = ^° (!^) 2 CW(A*/)Q(A*,) ( 8 ) 

where Q(Av) incorporates the effect of the Av dependence 
of O and(ff). We are mainly interested in the Av de- 
pendence, and we do not show the v dependence explicitly. 
Equation (8) relates the angular power spectrum of the fore- 
ground contribution to the visibility correlations which can 
be determined from our observations. 

The system noise makes a contribution 

Ar 2 (U 1 ,fi;U2,i>a) = <5u 1 ,u 2 <W2 V W 2 ) (9) 

which is non-zero only when a particular visibility is cor- 
related with itself. For a single polarization, the rms. noise 
in the real part (or equivalently the imaginary part) of a 



visibility is expected to be (Thompson, Moran & Swenson 
1986) 

a = ^ fcB Z^l_ (10) 
A eff VAvAt 

where T aya is the total system temperature, ks is the Boltz- 
mann constant, A e ff is the effective collecting area of each 
antenna, Av is the channel width and At is correlator inte- 
gration time. For the GMRT parameters 6 this is predicted to 
be g = 1.03 Jy for a single polarization. We have combined 
both polarizations, and so the variance in each visibility of 
the final data that we have analyzed is 2 a 2 . In eq. (9) the 
variance of the real and imaginary parts of the noise in a 
visibility contribute in quadrature and we have {N 2 ) = 4<r 2 . 



3.1 Estimating the visibility correlation. 

We use the estimator 

V 2 (U, Av) = V(U, Vi)V*(U + AU, Vi + Av) (11) 

where the bar denotes an average over the data under the 
assumptions 

(i) The U dependence is isotropic ie. V2 depends only on 
the magnitude U and not the direction of U 

(ii) The Av dependence is the same if the frequency ori- 
gin Vi is shifted to another channel Vj in the observation 
frequency band. 

(iii) Only visibilities V(U + AU, Vi + Av) at baselines U+ 
AU within a disk of radius | AU |< D < (tyOo) -1 centered 
at U are correlated with V(U, Vi), and V% is averaged over 
this disk. 

Note that the second assumption above implies that 
Vi{U, Av) gives an estimate of the average Av dependence 
across the entire frequency band. It also implies an average 
over positive and negative Av values. Besides this, the esti- 
mator is averaged over bins in U (Ui — U2, Vi — U3, ...). so 
that we have ViiJJi, Av) at a few values Ui corresponding to 
the average baseline of the bins. 
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The correlation of a visibility with itself introduces a 
noise contribution in the expectation value of this estima- 
tor. The noise contribution can be avoided (eg. Begum, 
Chengalur & Bhradwaj 2006) by excluding self- correlations 
ie. the visibility V(U, Vi) is correlated with every baseline 
V(U + AU, vt + Av) within a disk |AU| < D except itself. 
The expectation value of the estimator has a value 



(V 2 (U,Av)) = F 2 (U,Av). 



(12) 



which provides an unbiased estimate of the foregrounds. The 
system noise makes a contribution only to the uncertainty or 
the error in the estimator. The expectation value of the es- 
timator is real. The value of the estimator determined from 
an observation will, in general, have a real and an imagi- 
nary part. The real part contains the foreground informa- 
tion, whereas the imaginary part of the observed value of 
the estimator can be attributed to statistical fluctuations in 
the foregrounds and the noise. 

3.2 Error Estimates 

The expected uncertainty or statistical fluctuations in the 
real part of the estimator 

VW 2 ) 2 > = yf {(&-{&))*) (13) 

is the sum of two contributions 

{(AV,) 2 } = (AF 2 ) 2 + (AiV 2 ) 2 . (14) 

If we assume that the foregrounds are a Gaussian random 
field, the foreground contribution to the error is 



[AF 2 (K,A^)] 2 = ] i- 



F 2 (Ui,Q) + F2{Ui,Av) 



(15) 



where Ne is the number of independent estimates of 
F2(U, Av) that contribute to V 2 (Ui,Av). The baselines 
within a disk of radius ~ (tt6o)~ in uv space (Figure 1) are 
correlated, and all the baselines within such a disk provide 
only one independent estimate of the visibility correlation. 
For each U bin Ne is determined by counting the number 
of such regions with the uv coverage of our observations. 

The system noise contribution in any two visibilities are 
uncorrelated, and hence 



(16) 



where Np is the number of visibility pairs that contribute 
to the estimator V 2 {U,Av) for a particular U bin and Av 
separation. 

The error in the imaginary part of the estimator also is 
a sum of two contributions. The foreground contribution is 
somewhat different from eq. (15) and we have 



[AF 2 ([/,,Az,)] 2 = ^- 



F 2 (Ui,0)-F 2 {Ui,Av) 



(17) 



while the system noise contribution is the same a eq. (16). 



4 FOREGROUND MODEL PREDICTIONS 

We consider only the two most dominant foreground com- 
ponents namely extragalactic radio sources and the diffuse 



synchrotron radiation from our own Galaxy. The free-free 
emissions from our Galaxy and external galaxies is around 
1% of the total foreground contribution (Shaver et al. 1999), 
and we ignore this in our analysis. For each foreground com- 
ponent the MAPS can be modeled as 

where Vf = 130 MHz, and for each foreground component 
A, P and a are the amplitude, the power law index of 
the angular power spectrum and the mean spectral index 
respectively. The actual spectral index varies with line of 
sight across the sky and this causes the foreground con- 
tribution to decorrelate with increasing frequency separa- 
tion Av — \ V2 — vi | which is quantified through the fore- 
ground frequency decorrelation function Ii{v\v 2 ) (Zaldar- 
riaga, Furlanetto & Hernquist 2004) which has been mod- 
eled as 



h(vi , v 2 ) = exp 



-logi 



(19) 



The model parameters values that we have used are dis- 
cussed below and are given in Table 1. 

Resolved extragalactic radio sources (point sources) 
dominate the radio sky at 150 MHz. Di Matteo et al. (2002) 
have used the 6C survey (Hales, Baldwin & Warner 1988), 
and the 3CR survey and the 3 CRR catalogue ( Laing , Riley 
& Longair 1983) to estimate this contribution. The limiting 
flux density of these surveys was ~ 100 mjy and the ex- 
trapolation to fainter sources is rather uncertain. Di Matteo 
ct al. (2002) have fitted the differential source counts using 
a double power-law with the change in slope occurring at 
880 mjy. Since the brightest source in our image has a flux 
density below 880 mjy we use only the fit to the fainter part 



dN _ 4000 / S V 
dS ~ Jy-Sr' \ Uy ) 



(20) 



These sources make two distinct contributions to MAPS, the 
first being the Poisson noise arising from the discrete nature 
of these sources and the second arising from the clustering of 
the sources. Table 1 shows the respective parameters based 
on the estimates of Di Matteo et al. (2002) who assume that 
these sources are clustered like galaxies today or as Lyman- 
brcak galaxies (Giavalisco et al. 1998 ) at z ~ 3. Using these 
in eq. (8) to calculate the foreground contribution to the 
visibility correlation at 153 MHz for Av — 0, we have the 
Poisson term 



F 2 (U,0) = 7.6 



and the clustering term 

, 0-5 

F 2 (U,0) = 0.51 



Jy 2 , 



(—) 

\ 1000/ 



(21) 



(22) 



Here it is assumed that sources with flux greater than S c 
have been identified from continuum images and removed 
from the data. The brightest source in our initial image has 
S ~ 890 mjy and we use this value for S c when comparing 
model predictions with results from data I. For data R we 
have used S c = 8mJy as we have used this as the limiting 
value for our source subtraction (Section 2). 
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Table 1. Fiducial values of the parameters used for characterizing 
different foreground contributions 



Foregrounds A(mK 2 ) a £ 

Point source 1.2 x 10 4 (%^-) 1,25 2.07 1 
(Poisson part) 

Point source 6.1 X 10 3 (h^ 4 ) ' 2.07 1.1 2 
(clustered part) 

Galactic synchrotron 700 2.80 2.4 4 



estimator V2(U, Av) averages positive and negative Av val- 
ues. We use B = 1 to make an order of magnitude esti- 
mate. The expected change in F2(U, Av) is ~ 3 x 1CP 2 % 
for Av = 2.5 MHz. The key point here is that F 2 (U,Av) is 
predicted to change very slowly with Av, and the change is 
also very small. 



5 RESULTS AND DISCUSSION 

We have determined the observed value VziU, Av) of the 
visibility correlation estimator V 2 {U, Av) for data I and 
data R which are before and after source subtraction re- 
spectively. Baselines in the range 20 < U < 2 x 10 4 , and 
frequency channels 21 to 100 were used for the analysis. Vis- 
ibilities V(U + AU, v + Av) within the disk | AU < D = 5 
were correlated with V(\J,v). Here Av was restricted to 
| Av |< 2.5 MHz which corresponds to a separation of 40 
channels. Note that the correlation of a visibility with itself 
was not included. The value of D was chosen such that it 
is both less than (-7r#o) _1 = 8, and also large enough that 
a reasonable number of visibility pairs that contribute to 
the correlation. Figure 5 shows Vi(U, Av) as a function of U 
for Av = 0. Equivalently, we may also interpret this as the 
multi- frequency angular power spectrum Ci(Av) at Av — 0. 

For both the data-sets the real part of V2(U, 0) is found 
to be considerably larger than the imaginary part. This is 
consistent with the discussion of Section 3.1, and we expect 
the real part to provide an estimate of the foreground contri- 
bution V2(U, 0). The 1 — a error bars shown in the figure have 
been determined based on the error estimates discussed in 
Section 3.2. The uncertainty in F2(U, 0) is mainly due to the 
limited number of independent estimates, the system noise 
makes a smaller contribution. Though the results for data 
I over the range 200 < U << 2 x 10 4 looks like a power law 
V 2 (U,0) oc U~ a with a very small slope < a < 0.25, we 
do not find a fit with an acceptable value of \ 2 P er degree 
of freedom. 

The real part of V2(U, 0) falls to nearly one- fourth of its 
original value at most of the U bins when the directly de- 
tected sources are subtracted out. This indicates that a large 
part of the contribution to V2(U, 0) in data I is from these 
resolved sources, and we may interpret V2{U, 0) as arising 
primarily from these sources. Data R is expected to con- 
tain contributions from point sources below the detection 
limit of our image, diffuse sources, system noise, limitations 
in our imaging and source subtraction procedure and resid- 
ual RFI. We will assume for the moment that these effects 
can be ignored, but return to this issue later in this section. 

Figure 6 shows the observed V2{U, 0) plotted against the 
predictions of the foreground models discussed in Section 4. 
The brightest source in our image has flux 890 mjy. Based 
on this we use S c = 900 mjy for the point source contribu- 
tion to data I. The clustering of point sources dominates 
at baselines U < 150 (9 > 0.7°), while the Poisson fluctua- 
tions of the point sources dominates at larger baselines. The 
diffuse Galactic synchrotron radiation is much smaller than 
the point source contribution at all baselines. The errors 
in the model prediction are quite large and are mainly due 
to the Poisson fluctuations of the point sources. The model 
predictions are found to be consistent with the observed 
values of V2{U, 0) except at the smallest U value which cor- 



The uncertainty or error in the model prediction for 
these radio sources is also a sum of two parts. The error in 
the clustering part can be estimated using eq. (15). For the 
Poisson part the variance of F2 involves the fourth moment 
of the differential source count and we have 



[AF 2 (U,0)] 2 



Jy 



63.2 - 1.54 



Jy 



(23) 



The diffuse Galactic synchrotron radiation is believed to 
be produced by cosmic ray electrons propagating in the mag- 
netic field of the Galaxy (Ginzburg & Syrovatskii 1969) .This 
has an angular power spectrum that scales as Ci ~ l~ 2A 
(Tegmark et al. 2000), though this slope (0) is rather uncer- 
tain. The analysis of radio surveys at 408 MHz, 1.42 GHz, 
and 2.326 GHz (Haslam et al. 1982; Reich 1982; Reich & Re- 
ich 1988; Jonas, Baart, & Nicolson 1998) show the spectral 
index to be a ~ 2.8 which is in general agreement with re- 
sult of Platania et al. (1998). For the synchrotron radiation, 
in Table 1 we have adopted the parameters from Santos et 
al. (2005) which gives 

F 2 ([/,0)=4.2x KT 3 (JL) 24 Jy 2 . (24) 

We note that the amplitude of the synchrotron contribution 
is very sensitive to the spectral index whose value is quite 
uncertain. The value is in the range 2.5 < a < 3, and the 
amplitude increases by nearly an order of magnitude if a — 3 
instead of a = 2.8 as assumed here. 

The error for the synchrotron prediction can be calcu- 
lated using eq. 15. The total error in the model predictions 
is calculated by adding the variances from the different con- 
tributions. 

For the frequency separations of our interest (Av < 
2.5 MHz), for all the foreground components the (v2/vf) a 
term in equation (18) introduces a larger Av dependence in 
Ci (Av) as compared to the frequency decorrelation function 
I(vi, V2)- When calculating F2(U, Av) it is necessary to also 
incorporate Q(Av) (eq. 8) which has the Av dependence 
arising from 9o and (8B/dT) v . All of these predict a smooth 
Av dependence, and we may use a Taylor series expansion 



F 2 (U,Av) = F 2 (U,0) 



'-(7)' 



(25) 



where B is a constant of order unity. The Av/v term does 
not appear in eq. (25). This term cancels out because the 
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U 



u 



Figure 5. This shows the real (upper curve) and imaginary (lower curve) parts of the observed visibility correlation V^(£/, 0) as a function 
of U for the two data-sets indicated in the figure. As shown here, this may also be interpreted as C;(0) as a function of I. 



responds to an angular scale of ~ 1.8°. At these baselines 
the convolution with the primary beam pattern (eq. (6)) 
becomes important. We have not included this, and the ac- 
tual model predictions would possibly be somewhat smaller 
if this were included. As noted in Section 4., the amplitude 
of the synchrotron contribution is very sensitive to the value 
of the spectral index. The amplitude decreases by a factor 
of ~ 18 if a — 2.5 instead of the value a = 2.8 used here. 
This changes the total foreground contribution only at small 
baselines (U < 100) where the model then becomes consis- 
tent with our observations. 

The limiting flux for source subtraction is ~ 8 mjy, and 
hence we use S c ~ 10 mjy for data R. The model prediction 
is dominated by Galactic synchrotron radiation at U < 150, 
point source clustering in the range 150 < U < 2 x 10 3 
and point source Poisson fluctuations at U > 2 x 10 3 . The 
model predictions fall short of the observations at all base- 
lines except the smallest U value where it overshoots the 
observations. Since the model prediction for S c ~ lOmJy 
falls very much short of the observations, we also consider 
S c = 100 mjy where the dominant contribution is Galac- 
tic synchrotron at U < 60, point source clustering in the 
range 60 < U < 400 and point source Poisson fluctuations 
at U > 400. We find that the observations are slightly above 
the 1 — a error-bars at baselines U > 100, whereas they ex- 
ceed the model predictions at baselines U < 100. A point 
to note is that at the smallest baseline the prediction for 
the Galactic synchrotron radiation exceeds the observation. 
This may be a consequence of the possibility that the back- 
ground radiation is relatively low in the direction of our 
observation. Estimates from the Haslam et al. (1982) map 
at 408 MHZ show a relatively low brightness temperature of 
~ 30 K towards the direction of our observation. 

We quantify the Av dependence of Vi{U, Av) using 
k(U,Av) which is defined as 



k(U, Av) = 



V 2 (U,Av) 



(26) 



V 2 (U,0) 

We expect the visibilities V^U, v) and V(U, v + Av) to 
get decorrelated as Av is increased, and hence we expect 



< k(U, Av) |< 1. Figure 7 shows k(U,Av) for different 
values of U. The foreground models predict a smooth Av 
dependence for k(U, Av). The departure from k(U, Av) = 1 
is predicted to be less than 1% for Av < 2.5 MHz. The ob- 
served behavior of k(U, Av) is quite different from the model 
predictions. At the small baselines U < 1000 we find that 
k(U, Av) falls sharply within the first three channels. In the 
U = 47 bin k(U, Av) fluctuates at large Av whereas it re- 
mains roughly constant at U = 360. In both cases this value 
of k(U, Av) is smaller for data R as compared to data I. At 
U = 2200, for data I k(U,Av) falls gradually with increas- 
ing Av, and the visibilities are uncorrelated (k(U, Av) ~ 0) 
by Av ~ 2.5 MHz. Interestingly, for data R we find that 
n(U, Av) shows a sudden increase to n(U, Av) > 1 at very 
small Av (< 0.5 MHz), after which k(U,Av) falls and be- 
comes negative by Av ~< 2 MHz. It appears that in this U 
bin our source subtraction procedure has introduced excess 
correlations between the visibilities at small Av and intro- 
duces anti correlations at large Av. At U = 4200, for data 

1 the value of k(U, Av) oscillates with increasing Av. At 
large Av data R also shows a similar behavior except that 
the k(U, Av) values are smaller. The behavior of data R is 
quite different from that of data I at very small Av where 
there are two small oscillations that cross k(U, Av) = 1. 

The first point that emerges from our results is that the 
observed visibility correlations V^(J7, 0) is consistent with the 
predictions of the existing foreground models at all baselines 
except the smallest one which probes angular scales ~ 1°. 
The observations are in excess of the model prediction at 
the smallest baseline. The second point is that VziJJ, Av) 
shows considerable Av dependence, there being changes of 
order unity within Av — 2.5 MHz. This rapid change in the 
visibilities Va(U,v) across frequency channels is contrary to 
the foreground models which predict changes less than 1%. 

It is well appreciated that accurate subtraction of the 
foreground emission requires very exacting calibration. In 
contrast, we have followed fairly standard calibration proce- 
dures. As such it seems likely that the discrepancy between 
our observations and existing predictions is probably not 
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u u 

Figure 6. The thick solid line shows the real part of the observed visibility correlation V2(U, 0) as a function of U for the two data-sets 
indicated in the figure. As shown here, this may also be interpreted as C;(0) as a function of I. For data I the thin solid line shows the 
total model prediction for S c = 900 mjy. Also shown are the contributions from point source Poisson (dash-dot), point source clustering 
(dot) and Galactic synchrotron (dash-dot-dot-dot). For data R the thin solid line shows the total model predictions for S c = 100 mjy 
and and the long dashed line for 10 mjy. The dash-dot-dot-dot curve shows the Galactic synchrotron contribution. 




Az/ MHz 

Figure 7. This shows k(U, Au) as a function of Au for the different U values shown in the figure. The upper curve (at large A^) shows 
data I while the lower shows data R. 



genuine; indeed there are a several purely instrument related 
possibilities that may account for the discrepancies between 
our observational findings and existing models for the fore- 
ground emission. We take up first the issue of calibration 
error which will introduce phase and amplitude errors in the 
visibilities. The fact that the values of k(U, Au) are generally 
smaller for data R as compared to data I may be inter- 



preted as indicating that the visibilities V(U, v) are a com- 
bination of two parts, a correlated part which arises from for 
e.g. the effect of calibration errors on discrete sources, and 
another whose contribution to different channels is uncorre- 
cted. The "halos" that we see around the bright sources is a 
clear indication that calibration problems exist in our data. 
Phase errors which vary with channel would cause decorre- 
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lation of the visibilities across different frequencies. Further, 
one would expect that the phase errors increase with in- 
creasing baseline length, which is qualitatively consistent 
with what we see in Fig. 7. In contrast to the situation 
for k(U, Az/), the contribution from the source subtraction 
residuals to V2(U, 0) (Fig. 6) can be estimated to be small 
as follows. There are only ~ 100 imaging artifacts with ab- 
solute value of flux > 20 mjy (Data R, Figure 3) , while 
about 10,000 such sources would be needed to produce the 
observed visibility correlation of ~ 4Jy 2 (Data R, Figure 
6). 

The 2D Fourier relation between the sky brightness and 
the visibilities assumed in Section 3 is not strictly valid for 
GMRT's large field of view (#fwhm = 3.8°). In addition to 
u — v which are the components of the baseline in the plane 
normal to the direction of observation, it is also necessary 
to consider w the component along the observing direction. 
This is a possible source of error in our visibility correla- 
tion analysis. To asses the impact of the w term we have 
repeated the analysis using only a limited range of baselines 
for which w < 100. We find that limiting the maximum w 
value does not make any qualitative change in our results. 
The conclusions are unchanged even if we impose w < 50. 

Residual RFI is another possibility. The visibilities were 
clipped at 12 Jy (Section 2.) and this is expected to remove 
the strong RFI, but weak RFI contributions will persist in 
the data. The RFI electric fields at any two antennas is cor- 
related with a time delay r which depends on position of 
the RFI source relative to the antennas and the direction of 
observation. The RFI contribution behaves like the system 
noise if r is greater than r c the coherence time of the RFI 
signal. In this case the RFI effectively increases a the rms. 
fluctuations of the visibilities. This only changes the error 
estimates, and does not affect the expected visibility corre- 
lations. RFI sources for which r < r c are expected to affect 
the visibility correlations. This contribution will depend on 
the distribution of the time delays rs and the frequency 
spectrum of the RFI sources. The analysis of this is beyond 
the scope of this paper. Work is currently underway at the 
GMRT to implement more sophisticated real time as well as 
offline RFI mitigation schemes. Future observations will help 
assess the improvement that these schemes as well as better 
calibration procedures make on the problem of foreground 
subtraction. Polarization leakage is another important issue 
that we plan to take up in future work. 
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